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Abstract. We provide two algorithms for computing the volume of the 
convex polytope 0,:= {x gWI\Ax< b}, for A £ R™'<", b G K". Both 
algorithms have a 0(n™) computational complexity which makes them 
especially attractive for large n and relatively small m, when the other 
methods with 0(m") complexity fail. The methodology which differs 
from previous existing methods uses a Laplace transform technique that 
is well suited to the half-space representation of fi. 



1. Introduction 

In this paper, we are interested in the exact computation of the volume 
of the convex polytope := {x G | Ax < b}, for some given matrix 
^ G ]R"^><" and vector beW. 

Computing the volume of a convex polytope Q is difficult. Basically, 
methods for exact computation of this volume can be classified according 
to whether one has a half-space representation of as above, or a vertex 
representation, that is, when 0, is given by its list of vertices (triangulation 
methods), or when both descriptions are available. For instance, Lasserre 
's algorithm ^ requires a half-space description, whereas Delaunay 's tri- 
angulation (see e.g. 1^) or Von Hohenbalken's simplicial decomposition |ll] 
require the list of vertices. On the other hand, both Lawrence 's formula [|] 
and Cohen and Hickey 's triangulation method Q require the double (half- 
space and vertex) description of the polytope. For an updated review of the 
above methods and their computational complexity, the interested reader is 
referred to Biieler et al [^. In particular, improved versions of some of the 
above algorithms are also described in Q. The computational complexity 
is also discussed in Dyer and Frieze ||^]. In a different spirit, Barvinok 
approximates the volume by computing the integral of exp (c, x) over for 
a small c. The latter integral reduces to evaluate at each vertex v ^l, the 
integral of exp (c, x) over the smallest convex cone at v, which contains 
Q. Interestingly, the latter integrals are computed via a Fourier transform 
technique. 

In general, when Q has a half-space representation, the methods described 
in Biieler et al @ have a computational complexity that is exponential in 
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n, the dimension of the underlying affine space. While those methods work 
well for relatively small n and possibly large m, they become very time- 
consuming and even fail for large (or even not so large) n. This was the 
motivation for an alternative method that could work in the "dual" context 
of possibly large n and relatively small m. 

Here we suppose given a half-space representation of Q. The alternative 
method that we propose is conceptually very simple (as well as the compu- 
tations involved) and differs from previous existing methods. The idea is to 
consider the volume of 17 = {x > 0; Ax < 6} as a function g : M™— >M-|_ of the 
right-hand side b G for which we provide a simple explicit expression of 
its Laplace transform G : C"— >C in closed form. It then suffices to apply the 
inverse Laplace transform to G, which, in the present context, can be done 
efficiently by repeated applications of Cauchy 's Residue Theorem for the 
evaluation of one-dimensional complex integrals. We propose and describe 
two such algorithms. 

As already mentioned, the 0(n™) computational complexity of both al- 
gorithms makes the method especially attractive for large n and relatively 
small m, when the other methods with computational complexity 0{m"') 
would fail. This method can also be viewed as "dual" of the latter methods 
which work in the original space M" with the matrix A, as we instead work 
in the space M™" of "dual" variables associated with the constraints, and 
the cone {u > 0, A'u > 0} (via the Laplace transform), which explains the 
computational complexity 0{n"^) (in lieu of 0{m"')). 

2. Main result 

Let Cj := (1, 1, • • • ) be the unit vector of W for i > I. Let y G and 
A G ]^'"X'T- real-valued matrices such that the convex polyhedron 



is compact, that is, i^{y) is a convex polytope. The notation stands for 
the semi-closed interval [0, oo) C M. 

Now consider the function g : M^^M defined by 



Jniy) 

and let G : C" ^ C be its m-dimensional two-sided Laplace transform, that 
is. 



(2.1) 



Q{y) := {xeRl\Ax < y} 



(2.2) 





(2.3) 




We have the following result : 



Theorem 2.1. Let 0,{y) be the convex polytope in (2.1), functions g and G 
are defined as in ^2.^ and ( \2.3l ) respectively, and assume that x = is the 
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only solution of the system {x > 0, Ax < 0}. Then : 

I 1 3f?(A) > 0, 

^{A'X) > 0. 



(2.4) G(A) = 
Moreover, 

(2.5) g{y) 



1 



ci+ioo 



where the real constants c > satisfies A'c > 0. 
Proof. Apply the definition ( ^ ) of G, to obtain : 



e<^'2'>G(A)(iA 



G(A) 



y>Ax 



dx 



x>Q, Ax<y 



dy 



dy 



dx 



1 



nr=i A. 
1 



e-<^'^'->dx, 
1 



3?(A) > 



with 



5R(A) 
K(A'A) 



> 
> 



And (|2.5| ) is obtained by a direct application of the inverse Laplace trans- 
form. It remains to show that, indeed, the domain {5R(A) > 0, K(A'A) > 0} 
is nonempty. However, this fact follows from a special version of Farkas' 
lemma due to Carver (see e.g. Schrijver [l^, (33), p. 95]), which (adapted 
to the present context) states that {u > 0, A'u > 0} has an admissible so- 
lution u € M™^ if and only if (x, y) = is the only solution of the system 
{Ax + y = 0, x>0, y> 0}. In other words, x = is the only solution of 
{x > 0, < 0}. □ 

Remark 2.2. A necessary and sufficient condition for fl{y) to be compact is 
that there exists some u G such that A'u > Cn- This is a consequence of 
the well-known Farkas Lemma. 



As mentioned above, computing g{y) via ( p. 5] ) reduces to computing the 
Laplace inverse of G(A). In our case, this can be done quite efficiently even 
for large size problems. We first slightly modify our problem as follows : 

Suppose that we want to compute the volume of the convex polytope 
{x > 0; Ax < b} with 6 > 0, that is, we want to evaluate g{y) at the point 
y := h ^ W^. We may and shall assume, without loss of generality, that 
yi = I for every i = 1, . . . m. Otherwise, just divide by > 0. 

The problem is then to compute the value h{l) of the function h : M_|_^M 
given by 



(2.6) h{z) = g{emz) 



1 



(27ri) 



ci+ioo 
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where the real vector < c G M"* satisfies A'c > 0. Computing the complex 
integral ( |2.6|) can be done in two ways that are explored below. We do it 
directly iii|^ by integrating first with respect to (w.r.t.) Ai, then w.r.t. A2, 
etc., or indirectly in by first computing the one-dimensional Laplace 
transform H oi h and then computing the Laplace inverse of H. 



3. The direct method 

To better understand the methodology behind the direct method and for 
illustration purpose, consider the case of a convex polytope i7 with only two 
(m = 2) nontrivial constraints. 

3.1. The m = 2 non trivial constraints example. Let A G M^^" be 
such that X = is the only solution of {x > 0, Ax < 0}. Moreover, suppose 
that A' := [a\b] with a,b M". For ease of exposition, assume that 

• ajbj 7^ and aj 7^ bj for all j = 1,2, . . . n. 

• aj/bj / ak/bk for all j, k = 1,2, . . . n 

Observe that these assumptions are satisfied with probability one if we add 
to every coefficient ai,bi a perturbation e, randomly generated under a uni- 
form distribution on [0, e], with e very small. 
Then : 

G(A) - ' ' / > 



A1A2 n"=i(«iAi + 6JA2)' \ 5R(aAi + 6A2) > ' 

Next, fix ci and C2 > such that ajci + bjC2 > for every j = 1,2, . . . n, 
and compute the integral ( |2.6[ ) as follows. We first evaluate the integral 



(3-1) h = — — r~rTTT^Ai, 

using classical Cauchy 's residue technique. That is, we: (a) close the path 
of integration adding a semicircle V of radius R large enough, (b) evaluate 
the closed integral using Cauchy 's Residue Theorem ^ Theor. 2.2, p. 112], 
and (c) show that the integral along T converges to zero when R ^ 00. 

Now, since we are integrating w.r.t. Ai and we want to evaluate h(z) at 
z = 1, we must add the semicircle T on the left side of the integration path 
K(Ai) = ci because e^"**^ converges to zero when Ai —00. Therefore, we 
must consider only the poles of G(Ai, •) whose real part is strictly less than 



ci (with A2 being fixed). Then, the evaluation of (3J.) follows easily, and 



1 " -Q-(bj/aj)zX2 

^ n •=! b, + fr^^ 6, A^ Uk^,i-akb,/a, + b,) ' 

Recall that ^{—\2bj/aj) = —C2bj/aj < c\ for each j = 1,2, ...n, and 
G(Ai, •) has only poles of first order (with A2 being fixed). 
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Therefore, 

Hz) 



2 /•C2+ioo qZX2 



27ri .Ic2-ioo ^2 

C2+ioo 



h dX2 



2"^* Jc2~ioD ^2^^ Y[j=l 



1 /■C2+«00 

E — 



^{l~bj/aj)zX2 



^2 ^i^i nfc^i(^fc«i - akbj 



These integrals must be evaluated according to whether (1 — bj/aj)y is 
positive or negative. Thus, recalling that z > 0, each integral is equal to 

- its residue at the pole A2 = < C2 when 1 — bj/aj is positive, and 

- zero if 1 — bj/aj is negative because there is no pole on the right side of 

3f?(A2) = C2. 

That is. 



(3.2) 



Hz) = -7 



n! 



1 



ibj 



E 

bj/aj<l 



6,r 



(ijbjUk^jibkaj - akbj) 



Observe that the formula is not symmetrical in the parameters o, b. This is 
because we have chosen to integrate first w.r.t. Ai; and the set {j \ bj/aj < 1} 
is different from {j \cLj/bj > 1}, which would have been considered had we 
integrated first w.r.t. A2. In the latter case, we would have obtained 



(3.3) 



Hz) 



z 
nl 



1 



n-=i«i ^^j;^^,ci,b,Uk^jiakbj 

which is (|3.2| ) by interchanging a and b. Moreover, moving terms around, 
we get for free the following identity 

{aj - bj)"" _ 1 1 



(bj 



ai 



bkttj) 



(3.4) 



lajbjUk^jibkaj - akbj) lYj 



--ibj 



3.2. The direct method algorithm. The above methodology easily ex- 
tends to an arbitrary number m of non trivial constraints. One evaluates 
the integral of the right-hand side of (|2.6| ) by integrating first w.r.t. Ai, then 
w.r.t. A2, and so on. The resulting algorithm can be described with a tree 
of depth m + 1 {m + 1 "levels"). Let < c G M™ be such that A'c > 0. 

- Level is the root of the tree. 

- Level 1 is the integration w.r.t. Ai and consists of at most n + 1 nodes 
associated with the poles Ai := pj, j = 1, . . . n -|- 1, of the rational function 
Hi rijC^'-^)/^' seen as a function of Ai only. By the assumption on c, 
there is no pole pj on the line K(Ai) = ci. By Cauchy 's Residue Theorem, 
only the poles at the left side of the integration path 3?(Ai) = ci, say pj, 
j € Ii, are selected. 
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- Level 2 is the integration w.r.t. A2. After integration w.r.t. Ai, each of 
the poles pj, j G /i, generates a rational function of A2, A3, . . . , Am, which, 
seen as a function of A2 only, has at most n + 1 poles i = 1, . . . n + 1, 
j £ Ii. Thus level 2 has at most (n + 1)^ nodes associated with the poles 
pf{j). Assuming no pole pf{j) on the line K(A2) = C2, by Cauchy 's Residue 
Theorem, only the poles pf{j), € h, located on the correct side of the 
integration path 3?(A2) = C2 are selected. 

- Level k, k < m, consists of at most (n + 1)'^ nodes associated with 
the poles Pg{ii,i2, ■ ■ ■ ik-i), {h,i2, ■ ■ ■ ik-i) ^ h-i, s = 1, . . . n + 1, of some 
rational functions of A^, . . . , Am, seen as functions of A^ only. Assuming no 
pole on the line K(Afc) = c^, only the poles Pi^,{h, ■ ■ ■ ik-i), (ii, '^2, • • • ik) G 
Ik, located on the correct side of the integration path 3f?(Afc) = Ck, are 
selected. And so on. 

The last level m consists of at most (n + 1)™ nodes and the integration 
w.r.t. Xm is trivial as it amounts to evaluate integrals of the form 

(2^i)-i / AA-^^+i^e^^^-dAm, 



for some coefficients A, a, which yields A(az)'^ /nl for those a > 0. Summing 
up over all the nodes provides the desired value. 

Only simple elementary arithmetic operations are needed to compute the 
nodes at each level, as in Gauss elimination for solving linear systems. There- 
fore, the computational complexity is easily seen to be 0{n'^). 

However, some care must be taken with the choice of the integration paths 
as we assume that at each level k there is no pole on the integration path 
^{Xk) = Ck- This issue is discussed in §3^. The algorithm is illustrated on 
the following simple example with n = 2, m = 3. 

Example: Let $7 (262) C be the polytope 

0(2:62) := {x € \xi + X2 < z; —2xi + 2x2 < z ; 2xi — X2 < z}, 

whose area is 172^/48. 

Choose ci = 3, C2 = 2 and C3 = 1, so that ci > 2c2 — 2c3 and ci > C3 — 2c2. 

h{z) = j^J •../ e(^^+^^+^3).G'(A)<iA, 

with 



Ci— «oo J C^—lOO 



G{X) 



1 



AiA2A3(Ai - 2A2 + 2A3)(Ai + 2A2 - A3) ' 

Integrate first w.r.t. Ai; that is, evaluate the residues at the poles Ai = 0, 
Ai = 2A2 - 2A3 and Ai = A3 - 2A2 because < z, < ci, K(2A2 - 2A3) < ci 
and 3?(A3 — 2A2) < ci. We obtain 

hi^) = 7^=r-y2 / l2 + h + hdX2dX3, 



(27ri)2 



C3— 400 
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where 

_g{A2+A3)2 

" 2A2A3(A3-A2)(A3-2A2)' 

g(3A2-A3)2: 

" 6A2A3(A3-A2)(A3-4A2/3)' 

g(2A3-A2)2; 

" 3A2A3(A3-2A2)(A3-4A2/3)' 

Next, integrate I2 w.r.t. A3. We must consider the poles on the left side 
of K(A3) = 1, that is, the pole A3 = because 5R(A2) = 2. Thus, we get 
— e^'^2/4A2, and the next integration w.r.t. A2 yields — 

When integrating Is w.r.t. A3, we have to consider the poles A3 
A3 = 4A2/3, on the right side of 5?(A3) = 1; and we get 



A2 and 



-1 
Af 



322A2 2g2A25/3 



+ 



Recall that the path of integration has a negative orientation, so we have 
to consider the negative value of residues. The next integration w.r.t. A2 
yields z^{l- 25/48). 

Finally, when integrating 14 w.r.t. A3, we must consider only the pole 
A3 = 0, and we get e~^'^2/8A2; the next integration w.r.t. A2 yields zero. 
Hence, adding up the above three partial results, yields 



h{z) = 

which is the desired result. 



-1 



+ 1 



25 
48 



3.3. Paths of integration. In choosing the integration paths 3?(Afc) = c^, 
k = 1, . . . m, we must determine a vector < c € such that A'c > 0. 
However, this may not be enough when we want to evaluate the integral 
( p. 61 ) by repeated applications of Cauchy's Residue Theorem. Indeed, we 
have seen in the tree description of the algorithm (cf. § |3.2| ), that at each 
level A; > 1 of the tree (integration w.r.t A^), one assumes that there is no 
pole on the integration path 3?(Afc) = c^. 

For instance, had we set ci = C2 = C3 = 1 (instead of ci = 3, C2 = 2 and 
ci = 1) in the above example, we could not use Cauchy's Residue Theorem 
to integrate I2 or because we would have the pole A2 = A3 exactly on 
the path of integration (recall that 5?(A2) = ^?(A3) = 1); fortunately, this 
case is pathological as it happens with probability zero in a set of problems 
with randomly generated data A G i^»">^"- and, therefore, this issue could be 
ignored in practice. However, for the sake of mathematical rigor, in addition 
to the constraints c > and A'c > 0, the vector c G M*" must satisfy addi- 
tional constraints to avoid the above mentioned pathological problem. We 
next describe one way to proceed to ensure that c satisfies these additional 
constraints. 

7 



In §3.2 we have described the algorithm as a tree of depth m (level i being 
the integration w.r.t. Aj) where each node has at most n + 1 descendants 
(one descendant for each pole on the correct side of the integration path 
3^(Aj) = Cj). The volume is then the summation of all partial results obtained 
at each leaf of the tree (that is, each node of level m). We next describe how 
to "perturbate" on-line the initial vector c G M"* if at some level k of the 
algorithm there is a pole on the corresponding integration path 5ft(Afc) = c^. 

- Step 1. Integration w.r.t. Ai. Choose a real vector c := (c{, • • • , c^) > 
such that A'c > and integrate ( p.6| ) along the line 5R(Ai) = c\. From 
Cauchy 's Residue Theorem, this is done by selecting the (at most n+1) poles 
Ai := pj, j € /i, located on the left-side of the integration path 5R(Ai) = 
c\. Each pole Pj,j = l,...n + 1 (with pj := 0) is a linear combination 

/?g)A2 + . . . + /3]2a„ with real coefficients {/j]^^}, because A is a real- valued 
matrix. Observe that by the initial choice of c, 

m 

6, := min |ci - ^ > 0. 

k=2 

- Step 2. Integration w.r.t. A2. For each of the poles Pj, j & h, selected 



at step 1, and after integration w.r.t. Ai, we now have to consider a rational 



function of A2 with at most n + 1 poles A2 := pj{j) '■= Y1T=3 f^ik U)'^^ 



i = l,...n+ 1. If 



62 := min min 14 - V /J^^jlcfcl > 0, 

7G/1 1=1, ...n+1 ^ — ' 
k=3 

then integrate w.r.t. A2 on the line 5ft(A2) = c^. Otherwise, if ^2 = we set 
C2 := C2 + €2 and c| := for all A; 7^ 2, by choosing £2 > small enough to 
ensure that 

(a) A'c^ > 

m 

(b) 52 := min min |ci - V/3jf^(j)c^.| > 

jG/i 2=1, ...n+1 ^ 

(c) max |/3,-2^e2| < Si 

j=l,...n+l -' 

The condition (a) is basic whereas (b) ensures that there is no pole on 
the integration path 5ft(A2) = c^. Moreover, what has been done in step 1 
remains valid because from (c), cf — Y1T=2 f^jk '^t same sign as c} — 

Yjk=2 f^jk'^k^ ^^'^^ therefore, none of the poles pj has crossed the integration 
path 5ft(Ai) = c\ = cf, that is, the set Ii is unchanged. 

Then integrate w.r.t. A2 on the line 9?(A2) = c|, which is done via Cauchy 
's Residue Theorem by selecting the (at most (n+1)^) poles pfij), {j, i) € I2, 
located at the left or the right of the line K(A2) = c^, depending on the sign 
of the coefficient of the argument in the exponential. 



- Step 3. Integration w.r.t. A3. Likewise, for each of the poles pf{j), 
€ I2, selected at step 2, we now have to consider a rational function of 
A3 with at most n + 1 poles p^{j, i) := Y1T=4 /^Ife^i' O^fc, s = 1, • • • n + 1. If 



min min - ^ pfj {j, i)cl \ > 0, 

0,«)e/2 s=l,...n+l ^ 



then integrate w.r.t. A3 on the line 5ft(A3) = C3. Otherwise, if 5^ = 0, set 

C3 := C3 + £3 and cf. c\ for all A; 7^ 3, by choosing £3 > small enough to 
ensure that 

(a) A'c- > 

m 

(b) <53 := min min \cl-^p'f^{j,i)4\ > 

0,2)6/2 «=lvri.+l ^ 

(2) 

(c) max max |/?)3^(i)e3l < 62 

j£l\ 1=1, ...n+l 

(d) . max |/?g)e2 + /3j3^e3| < '^i 

]=l,...n+l ■' ■' 

As in previous steps, condition (a) is basic. The condition (b) ensures that 
there is no pole on the integration path 5R(A3) = C3. Condition (c) (resp. 
(d)) ensures that none of the poles (j) considered at step 2 (resp. none of 
the poles pj considered at step 1) has crossed the line 5R(A2) = = (resp. 
the line 5ft(Ai) = cf = c|). That is, both sets 1\ and I2 are unchanged. 

Then integrate w.r.t. A3 on the line 3f?(A3) = C3, which is done by selecting 
the (at most (n + 1)^) poles ps{j,i), {j,i,s) G 13, located at the left or 
right of the line 5R(A3) = C3, depending on the sign of the argument in the 
exponential. 

And so on. It is important to notice that the e^'s and c^'s play no (nu- 
merical) role in the integration itself. They are only used to (i) ensure the 
absence of a pole on the integration path 9f?(Afc) = c^, and (ii) to locate the 
poles on the left or the right of the integration path. Their numerical value 
(which can be very small) has no influence on the computation. 

4. The associated transform algorithm 

An alternative to the direct method permits to avoid evaluating integrals 
of exponential functions in ( |2.6| ) by making the following simple change of 
variable. Let Ai = p — Yl%=2 ^^^^ ~ ^j=i '^i ™ (|2.6|), so that 

rd+ioo 



^^'^ = (^ 

where 



C2+«00 



C2~lOO 



e^PG dp 

d—ioo 



d\2 ■ ■ ■ dXm, 



(4.1) G = G(p-^A„A2,...,A^). 
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We can rewrite h{z) as 



(4.2) h{z) = / e'PH{p)dp, with 

27ri Jd-ioo 

(4.3) H{p) := j^-^ •••/ GdA2...dA^.. 

Recall that G(A) is well defined on the domain K(A) > and ^{A'X) > 0; 
moreover, the real vector c is taken in this domain. Hence, the domain of 
definition of H(p) is given by the condition 

m 

mp)-Y,Cj,C2,...,Cm) G {y e I y > 0, A'y > 0}. 

i=2 

On other hand, recall that the system {x > 0, Ax < 0} has only one 



solution X = (see the hypotheses of Theorem 2.1). Hence, the function 
h{z) is identically equal zero when z < (see (2^) and (^.6|)). Therefore, 
H{p) is the one-sided Laplace transform of h{z). Moreover, it is also easy to 
see that there exists a real constant C such that h{z) = z"C/n! when z > 0. 
Therefore, 

H{p) = C/p''+^ 

and the main problem completely reduces to evaluating the constant C = 



h{l)n\ by integrating G in (4.3). 

Notice that we only need to evaluate m — 1 integrals. The function H{p) 
is called the associated transform of G(A). 

Again, the integral (|4.3| ) can be computed via repeated applications of 
Cauchy 's Residue Theorem (and as in the direct method algorithm of §^, 
some care is needed with the domain of integration and the location of the 
poles of G). The method is illustrated on the same example of two non 
trivial constraints (m = 2) already considered at the beginning of §^. 

4.1. The m = 2 non trivial constraints example. Let A G M^^" such 
that X = is the only solution of {x > 0, Ax < 0}. Write A' := [a\b] with 
a,b & M". To compare with the direct method, and as in the beginning of 
assume that ajbj 7^ for all j = 1, . . . and aj/bj ^ a^/bk for all j ^ k. 
Then : 

1 1 5R(A)>0, 

^ ^ A1A2 U]=MXi + bjX2y 5R(^'A) > 0. 

Fix A2 = p — Xi and choose a real constant ci > such that the system of 
inequalities 5ft(p) > ci and (aj — bj)ci + bj^{p) > for all j = 1, . . . n has a 
solution. We already know that there is at least one vector u = (ci, 5?(p) — ci) 
such that u > and A'u > 0. We obtain H{p) by integrating G{Xi,p — Ai) 
w.r.t. Ai, which yields 
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ci+ioo 



1 



1 



Next, we need to determine which poles of G{Xi,p — Ai) are on the left 
(right) side of the integration path 3?(Ai) = ci in order to apply Cauchy's 
Residue theorem. Let J+ = {j\aj > bj}, Jq = {j\aj = bj} and J_ = 
{j\aj < bj}. Then, the poles on the left side of 9f?(Ai) = ci are Ai = and 
Ai = —bjp/{aj — bj) for all j S J+ because —bj^{p)/{aj — bj 
the poles on the right side of 9?(Ai) = ci are Xi = p and Ai = 
for all j G J-. Finally, notice that G{Xi,p — Xi) has only poles or first order. 

Hence, computing the residues of poles on the left side of K(Ai) = ci, 
yields 



< Cj. Besides, 



-bjp/{aj 



H{p) 



1 



1 



+ 



-{aj - bj)''-\-^o\ 



+ Yl^ p2ajbj Uk^j,,^k^j{-pbjak + pajbk) 



p 



,n+l 



(^jbj Ilkj^ji^^jbk - bjak) 



aj/bj>l 



and one retrieves (3^) when we take Jo to be an empty set, in other words, 
when its cardinality | Jo| =0. Now, computing the negative value of residues 
of poles on the right side of 3?(Ai) = ci (we need to take the negative value 
because the path of integration has a negative orientation), yields 



H{p) 



p 



n 



i=i «i 



{b, - a,r 



■^^^ ajbj Uk^jibjak - ajbk) 



bj/aj> 



and we also retrieve 



4.2. The associated transform algorithm. As for the direct method al- 
gorithm, the above methodology easily extends to an arbitrary number m 
of nontrivial constraints. The algorithm also consists of m (one-dimensional 
integration) steps. At each step, the several one-dimensional complex inte- 
grals are evaluated by application of Cauchy 's Residue Theorem [^, Theor. 
2.2, p. 112]. For same reasons as in the direct method, the computational 
complexity is easily seen to be 0{n^). 

The general case is better illustrated on the same example as in |3.2| . 
Again, to avoid the case of poles on the integration path in pathological 
examples, some care is needed when one specifies the integration path at 
each step of the algorithm. 



Let ^1(2:62) 
0(2:62) :-■ 



C M be the poly tope 

= {x € \ x\ + X2 < z; —2xi + 2x2 < z ; 2xi — X2 < z}, 
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whose area is 172:^/48. 

We can choose A3 = p — A2 — Ai and ci = C2 = 1 such that > 2, 

23?(p) > 5 and ^{p) < 5; and so 



1 rl+ioo /"1+ioo 

H^P) = Tr^ / M(A,p)dAidA2, 

\^l-Kl) Ji_joo Jl-ioo 



with 

M{\p) 



A1A2 (p - Ai - A2)(2p - Ai - 4A2)(2Ai + 3A2 - p) ' 



We first integrate w.r.t. Ai. Only the real parts of the poles Ai = and 
\i = {p — 3A2)/2 are less than 1. Therefore, the residue of the 0-pole yields: 



1 



l+ioo 



^^■^^ 2T^ij^.,oo A2(p-A2)(2p-4A2)(3A2-p)'^^'' 

whereas the residue of the {p — 3A2)/2)-pole yields 

1 j-l+ioo ^ 

^^■^^ 2^A-ioo A2(p-3A2)(p + A2)(3p-5A2)'^^'- 

At this point, we have to be careful; observe that 5/2 < < 5. How- 
ever, we cannot put = 3 because otherwise we will have a pole in the 
path of integration of (|4.4|) and (4^). We thus fix 3 < "^{p) < 5. Applying 
again Cauchy's Residue Theorem to (^3) at the pole A2 = (the only one 
whose real part is less than one), yields — l/2p^. 

Similarly, applying Cauchy 's Residue Theorem to ( [4. 51) at the poles A2 = 
and A2 = —p (the only ones whose real part is less than one), yields 
4/3p3 - l/8p3. 

We finally have that H{p) = (4/3 - 1/8 - l/2)/p^ = 17/24p3, and so 
h{z) = 17z^/48, the desired result. 

Concerning the pathological case of some poles on the integration paths 
at some step of the algorithm, the same remarks and similar remedies as for 
the direct method are valid (cf. §p.3[). 



5. Conclusion 

We have presented two algorithms for the exact computation of the vol- 
ume of a convex polytope given by its half-space representation. The method- 
ology behind both algorithms is conceptually simple as it reduces to invert 
the Laplace transform of the volume (considered as a function of the right- 
hand-side). Both algorithms are relatively easy to implement (with special 
care for the choice of the integration paths of the repeated one-dimensional 
integrals). Their 0(n™) computational complexity can make them espe- 
cially attractive for large n and small m, when the other methods (with 
half-space representation of $7) fail because of their 0{m") computational 
complexity. 
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